nN scattering and electromagnetic corrections 
in the perturbative chiral quark model 



V.E. Lyubovitskij a , Th. Gutsche a , Amand Faessler a , 

R. Vinh Mau b 

q \ & Institut fur Theoretische Physik, Universitdt Tubingen, Auf der Morgenstelle 14, 

CN ' D-72076 Tubingen, Germany 

P-ij h Laboratoire de Physique Theorique des Particules Elementaires, Universite P. et 

' M. Curie, 4 Place Jussieu, 75252 Paris Cedex 05, France 

(N 
> ■ 

■^J- , Abstract 

m 

We apply the perturbative chiral quark model to give predictions for the electro- 
magnetic 0(p 2 ) low-energy couplings of the ChPT effective Lagrangian that define 
the electromagnetic mass shifts of nucleons and first-order (e 2 ) radiative corrections 
to the nN scattering amplitude. We estimate the leading isospin-breaking correc- 
tion to the strong energy shift of the ir~p atom in the Is state, which is relevant for 
the experiment "Pionic Hydrogen" at PSI. 
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M. 

In [1,2] Weinberg and Tomozawa derived a model-independent expression for 
the S -wave ttN scattering lengths using the current algebra relations and the 
PCAC assumption. To reproduce the result for the nN scattering lengths one 
can use the specific Lagrangian with the nucleon field N refered to as the 
Weinberg- Tomozawa (WT) term [3]- [7] which is part of the effective Wein- 
berg Lagrangian. The effective Weinberg Lagrangian can be derived from the 
original cr-model [8] by performing a chiral-neld dependent rotation on the 
nucleon field [3]. On the quark level the same exercise was done in the frame- 
work of the cloudy bag model [9,10]. The chiral transformation eliminates 
the nonderivative coupling of the chiral (pion) field with the nucleons/quarks 
and replaces it by a nonlinear derivative coupling (axial vector term + WT 
term + higher order terms in the chiral field). Note, that both realizations of 
chirally-symmetric Lagrangians (the original cr-model and the Weinberg type 
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Lagrangian) should a priori give the same result for the nN S-wave scatter- 
ing lengths. In ref. [11] in the framework of perturbative chiral quark model 
(PCQM) [12,13] we demonstrate that the equivalence between the two theo- 
ries with nonderivative and derivative coupling of the chiral field to the quarks 
is also valid when including the photon field. 

The purpose of this letter is to calculate first-order (e 2 ) radiative correc- 
tions to the nucleon mass and the pion-nucleon amplitude at threshold. We 
thereby predict the 0(p 2 ) electromagnetic (e.m.) low-energy couplings (LECs) 
originally defined in the effective Lagrangian of Chiral Perturbation Theory 
(ChPT) [5,6]. Quantitative information about these constants is important 
for the ongoing experimental and theoretical analysis of decay properties of 
the n~p atom (for a detailed discussion see Ref. [14]). In particular, we give 
a prediction for the leading isospin-breaking correction to the strong energy 
shift of the n~p atom in the Is state. 

Following considerations are based on the perturbative chiral quark model 
(PCQM), a relativistic quark model suggested in [12] and extended in [13] 
for the study of low-energy properties of baryons. The model includes rela- 
tivistic quark wave functions and confinement as well as the chiral symmetry 
requirements. The quarks move in a self-consistent field, represented by scalar 
S(r) and vector V(r) components of a static potential with r = \x\ providing 
confinement. The interaction of quarks with Goldstone bosons is introduced 
on the basis of the nonlinear a-model [8]. The PCQM is based on the effective, 
chirally invariant Lagrangian C inv [13] 

C inv {x)=i){x){i 0- 7 V(r)- S(r) 
F 2 

+ — Tr[d fl Ud' x U^} 

where is the quark field, U = exp[i&/F] is the chiral field and F = 88 
MeV is the pion decay constant in the chiral limit [4,13]. In the following we 
restrict to the SU (2) flavor case, that is $ — > ii = nr. For small fluctuations 
of the mesons fields one can use the perturbation expansion in powers of the 
parameter 1/F. The PCQM was successfully applied to a-term physics and 
extended to the study of electromagnetic properties of the nucleon [13]. Similar 
perturbative quark models have also been studied in Refs. [15]. 

The quark field ip we expand in the basis of potential eigenstates as 

^i x ) = J2 b aU a (x) exp(-iS a t) - J2 dp v p(x) exp(i£ p t) (2) 
where the sets of quark {u a } and antiquark {vp} wave functions in orbits a and 
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(5 are solutions of the Dirac equation with the static potential. The expansion 
coefficients b a and dt are the corresponding single quark annihilation and 
antiquark creation operators. 

The direct way to generate the WT term in the Lagrangian (1) is through 
introduction of a unitary transformation on the quark field ip. The technique 
was, for example, performed in the context of the cloudy bag model [9]. With 
the unitary chiral rotation if) — > exp{— i'y 5 ^/ (2F)}i(j the Lagrangian (1) trans- 
forms into a Weinberg-type form C w containing the axial-vector coupling and 
the WT term: 

C w {x) =C (x) + CY ;s '» + o(tt 2 ), (3) 
Cq(x) = - S(r) - -y°V(r) - X - n(x) (□ + M 2 ) n(x), 

CY'' str (x) = 7^d„n(x) ^{x) Y 7 5 t^{x) - g§ d^{x) $(x) ^ r k ^(x), 

where CY ' str is the 0(tt 2 ) strong interaction Lagrangian, □ = <9 M <9 M and M v is 
the pion mass. 

In ref. [11] we demonstrate explicitly for the 7rN amplitude up to order (1/F 2 ) 
that the two effective theories, the original one involving the pseudoscalar cou- 
pling and the Weinberg type, are formally equivalent, both on the level of the 
Lagrangians and for the matrix elements. This equivalence is based on the uni- 
tary transformation of the quark fields, where, in addition, the quarks remain 
on their energy shell. The same relation also holds in a fully covariant for- 
malism, when quarks/baryons are on their mass shell. Particularly, we show 
that the Weinberg- Tomozawa result can be reproduced with the use of the 
original Lagrangian (1) if: i) we use the expansion of the chiral field up to 
quadratic terms and ii) we employ the full quark propagator including the 
antiquark components. The two forms of the Lagrangian also yield the same 
results when including the photon field. For the equivalence to hold it is essen- 
tial that the photons are introduced consistently in both formalisms, that is 
by minimal substitution. One can prove that both Lagrangians yield the same 
results for radiative corrections to the irN scattering amplitude at threshold. 

In this letter we apply the developed formalism to study e.m. corrections of 
nucleon properties, such as the mass and the nN scattering amplitude. We 
perform all calculations using the technically more convenient Lagrangian (3). 
Introduction of the e.m. field is accomplished by minimal substitution into 
Eq. (3): 

-> D^ip = + ieQA^ip, -> = + ee^A^j (4) 

where Q is the quark charge matrix. 
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Following the Gell-Mann and Low theorem [16] the e.m. mass shift Arrfff 1 of 
the nucleon with respect to the three-quark ground state \<f> > N is 



Am™ = N (M - \\ S(x°)d 4 x J dST[C em {x)C em {y)} |0 o >f (5) 

to order e 2 in the e.m. interaction. Subscript "c" in Eq. (5) refers to contribu- 
tions from connected graphs only. Superscript "iV" indicates that the matrix 
elements have to be projected onto the respective nucleon states. These nu- 
cleon states are conventionally set up by the product of single quark SU(6) 
spin-flavor and SU(3) C color w.f. (see details in [13]), where the nonrelativistic 
single quark spin wave function is replaced by the relativistic ground state 
solution. With the quark-photon interaction defined by the Lagrangian 



C em {x) = - e A^(x)Q^(x), (6) 

the e.m. mass shift Arrf^ 1 is generated by two diagrams: one-body (Fig. la) 
and two-body diagram (Fig. lb). 

The leading e.m. corrections (up to order e 2 /F 2 ) to the irN scattering ampli- 
tude at threshold are generated by the interaction Lagrangian 



Cf{x) = CY ]Str \x) + CY' em (x) (7) 

where £^' str i s given in Eq. (3) and the additional e.m. part £j^' em i s given 
by 



CY' em (x)=C em (x) + ^A^x)^(x)r[7r\x)T 3 -7T(x)rn%x)}r { .r) (S) 
-eA^{x)e U j 

The irN amplitude in the presence of 0(e 2 ) radiative corrections is given by 



ir^d^ix) - ^p-^(x)Yl 5 rMx) 



N (4>0; TTil £ ~ } I d 4 X! ... f d A X n T[CY( Xl ) . . . Cf{x n ) ] |0 O ; 7T,)f . (9) 
n=l n - J J 

The diagrams for 0(e 2 /F 2 ) radiative corrections to the nN amplitude at 
threshold are shown in Fig.2. To evaluate the diagrams in Figs.l and 2 we 
use the photon propagator in the Coulomb gaugeQ to separate the con- 
tributions from Coulomb and transverse photons. 

1 It can be shown that the results do not depend on the choice of the gauge. 
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First, we analyze the e.m. mass shift of the nucleon. The contributions of 
diagrams Fig. la and Fig. lb are given by 



Amr = e 2 • N (<j> \ J d 4 x J d 4 y 5(x°) D^(x - y) (10) 

x $ {x) r Q iG^x, y) Y Q My) \<f>o) N 
Amr 6 = y • N {M / d*x J d'ySix^D^x-y) 

x Mx) Y q Mx) My) i v Q My) \<i>o) N . 

where iGf(x,y) =< 0\T{ip(x)ip(y)}\0 > is the quark propagator in a binding 
potential. In the following we truncate the expansion of the quark propagator 
to the ground state eigen mode: 



iG4x,y)^iG (x,y) = u (x)u ($e- i£ ^°-^e(x -y ), (11) 

that is we restrict the intermediate baryon states to iV and A configurations. 
Inclusion of excited baryon states will be subject of future investigations. With 
the use of approximation (11) Am^™'" and Am^J™' 6 reduce to 



^T a = ^(N\jjQ^\N) J ||{^(-g 2 )] 2 -^[G p M (-g 2 )] 2 } 
Am %* = ^s / Cf {(N\j:Q {i) Q U) \N)[G P E (-q 2 )} 2 (12) 

3 — » 2 

-(N\J2Q^a^\N)^[GU-f)} 2 }, 

where \N) is the SU(6) spin-flavor w.f. of the nucleon. Here we introduce the 
proton charge (G P E ) and magnetic (G P M ) form factors (f.f.) calculated at zeroth 
order [13] (meson cloud corrections are not taken into account) with 



X ] n ,XN s G E {-f) = N {4> \ / d 3 xM?)l°M?)e i n<Po) N , (13) 



A 2m N XNfiM{ ~ f) = I d ^Mx)^M^)e i n<Po) N 

where xn s is the nucleon spin w.f. and &n is the nucleon spin operator. Note 
that the contributions of Coulomb and transverse photons to the e.m. mass 
shifts (see Eqs. (12)) are related to the nucleon charge and magnetic f.f., 
respectively. The sum 
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(N\J2(Q 2 ) {{) \N) + {N\J2Q {i:) Q U) \N) 



1=1 



1 for N = p 
for N = n 



(14) 



is equivalent to the charge matrix of nucleons (Qn being the nucleon charge). 
In the limit mjy — > oo (when we neglect the contribution of G P M in Eqs. (12)) 
we obtain for the e.m. mass shifts 



AmT = Am7 ;a + Am™* = ^ J |f [G|(-g 2 )] 2 (15) 

consistent with the result (Eq. (12.4)) of Ref. [17]. Hence, the e.m. mass shift 
of the neutron vanishes in the heavy nucleon limit. 

In the numerical analysis we use the variational Gaussian ansatz [13] for the 
quark ground state wave function with the following analytical form: 



uq(x) = N exp 



x 



2R 2 



ipax/R 



Xs Xf Xc 



(16) 



where N = [7r 3//2 i? 3 (l + 3p 2 /2)] -1 / 2 is a constant fixed by the normalization 
condition / d z xu\{x) uq{x) = 1; Xs, Xf-, Xc are the spin, flavor and color quark 
wave functions, respectively. Our Gaussian ansatz contains two model param- 
eters: the dimensional parameter R and the dimensionless parameter p. The 
parameter p can be related to the axial coupling constant g& calculated in 
zeroth-order (or the three quark-core) approximation: 



Therefore, p can be replaced by the axial charge gA by means of the matching 
condition (17). The parameter R can be physically understood as the mean 
radius of the three-quark core and is related to the charge radius of the proton 
in the leading-order approximation as 



(4>fo = / d*xul(x)x*u (x) = ^f^^jjr (18) 

In our calculations we use the value 5^=1.25 obtained in ChPT [4]. Therefore, 
we have only one free parameter, that is R. In the numerical studies [13] R is 
varied in the region from 0.55 fm to 0.65 fm, which corresponds to a change 
of {r 2 E )^ from 0.5 to 0.7 fm 2 . The exact Gaussian ansatz (16) restricts the 
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potentials S(r) and V(r) to a form proportional to r 2 . They are expressed in 
terms of the parameters R and p (for details see Ref. [13]). 



Using (16) the proton f.f. at zeroth order are determined as [13]: 



Ge{-cl 2 ) = exp 



q 2 R 2 



q 2 R 2 



1 - 



-K 



G M (-q 2 ) = exp^-?-^—^2m N R ]J k(i - -nj, 
With Eq. (19) the e.m. mass shift is finally given as 



(19) 



K= 2~To 9A - 



Am™ = —^= 

Am e n m = ^_ . 



i k 3 2 34 , , 

1 1 K K 1 K 

2 16 

-k[ 1 K 



(20) 



where a = 1/137 is the fine structure coupling. For our set of parame- 



ters g\ = 1-25 and R = 0.6 ± 0.05 fm we get Am: 



em 
P 



0.54 ± 0.04 MeV, 



Am em = _ 26 ± o.02 MeV and Am e n rn - Am™ = -0.8 ± 0.06 MeV. These 
and the following uncertainties in our results correspond to the variation of 
the parameter R. Our predictions are in qualitative agreement with the re- 
sults obtained by Gasser and Leutwyler using the Cottingham formula [17]: 
Am™ = 0.63 MeV, Am^ = -0.13 MeV, Am™ - Am™ = -0.76 MeV. 
To compare our prediction for the e.m. mass shifts of the nucleons with the 
result of ChPT [6], we recall the part of the ChPT Lagrangian [6] which is 
responsible for radiative corrections 



C' ChpT = e'N /, 1 



,,,,,, P-W . hi. ^-^r- t 



F 2 



2F 2 



The 0{p 2 ) low-energy constants (LECs) /i, f 2 and / 3 contain the effect of the 
direct quark-photon interaction. Matching our results for the nucleon mass 
shifts to the predictions of ChPT [6] with 



Am e p m \chPT = -±™U + h + v)> A < m \chPT = -AnaU + f 3 - (22) 



V 2/' " 1 V 2 

we obtain following relations for the coupling constants f\, ji and fz- 



h 



2 j R(2tt) 3 / 2 



29 



89 



18 



48 



(23) 



7 



fi + fs 



4R(2tt) 3 / 2 



125 473 o 



Our numerical result for f 2 = —8.7 ± 0.7 MeV is in good agreement with the 
value of f 2 = —8.3 ± 3.3 MeV [6,14] extracted from the analysis of the elastic 
electron scattering cross section using the Cottingham formula [17]. For /1 + /3 
we get -1.5 ±0.1 MeV. 

We furthermore give a prediction for the separate values of f±, f 3 and the 
ratio /1//2 as deduced from our model analysis of e 2 corrections to the nN 
amplitude. We denote the corresponding matrix element associated with the 
nucleon flavor transition Ni — > N 2 by M^J^ 1 . In the Coulomb gauge only six 
diagrams (Fig.2a-2f) contribute to the radiative correction to the nN ampli- 
tude at threshold. The contribution of the other diagrams (Fig.2g-2o) vanishes. 
The contributions of the different diagrams of Fig. 2 are as follow: 



M 



mN-2 



N 1 



a+b 



F 2 



v.; / d A x J d^yD^ix - y)ip (x)Y 
x (T^G4x,y)Q + QG4x,y)T>)YMy)\<f>o) N 
for Fig.2a and 2b where T ij = 25 ij r 3 - 5 i3 T j - 5 j3 r\ 



(24) 



M 



(e 2 );ij 

NiN 2 



for Fig. 2c, 



re 
F 2 



N , 



(0o I d x d yD^x - y)ip (x)Y 



xr^ ^)My)YQMy)\<Po) 



N 



(25) 



M 



(e 2 );ij 



N 1 



/r2 J d^x J d^yD^x -y)i) Q (x)Yrf 



(26) 



for Fig. 2d and 2e, 



M 



mN-2 



^ • N (<f>o\ / d 4 x J d 4 yD^(x-y)Mx)n 5 



(27) 



for Fig.2f. 



Truncating the quark propagator to the ground state mode the nN scattering 
amplitude at threshold including first-order radiative corrections is 
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Ml 



e 2 nN 



+ 



(An) 3 
114 dKq 2 ) 

1 1 

~8#(2tt) 3 / 2 



19q 



q I 



nN 



6m 2 N 



\G P E {-f)f 



(28) 



5<f 



\G P M {-f)f 



18m% 



953 



41 115 



K 



N 



29 89 2 
1 -l8 K + S K 



where 



Ana 



F 2 



N{ti 2 - (7i ) 2 }N and MJ 



wN 
2 



Aira 



N{n 2 r 3 - (7?f)7r°}iV 



and 



d±(q' 



q 2 R 2 



1 ± 



K 



A 1 -2k 

The contribution of the Coulomb photons to the amplitude Mf^ N is parame- 
trized by the proton charge form factor (Ge), transverse photons are related 
to the proton magnetic {Gm) and axial nucleon (Ga) f-f- where the latter is 
given by [13] 



G A {-q 2 ) = g A e W [- q ^)d.{q 2 ). (29) 

Again, as in the case of e.m. mass shifts, the amplitude M^ N is gauge- 
independent. In ChPT the corresponding amplitude is given by [6] 



Mt: N \cKPT = fiMf + ^Mf. (30) 

Comparing of Eqs. (28) and (30) we get the same expression for f 2 as already 
obtained from the e.m. mass shift (23). We also deduce the following relations: 



fi 

h 



8 j R(2tt) 3 / 2 
1 



41 

3 



115 

— K + 



953 



m(2ixf/ 2 



35 

y 



785 

/ 

18 



16 
1913 
"48~ 



-K 



-K 



The predicted ratio for /1//2 depends on only one model parameter p (or k) 
which is related to the axial nucleon charge g A calculated at zeroth order. 
In addition, the constants j\, ji and f'3 depend on the size parameter R 
of the bound quark. For our "canonical" set of parameters, g A = 1-25 and 
R = 0.6 ± 0.05 fm, used in the calculations of nucleon e.m. form factors and 
meson-baryon sigma terms [13] we obtain: 
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fa = -19.5 ± 1.6 MeV, 



fa = -8.7 ±0.7 MeV, 



(31) 



fa = 18 ±1.5 MeV, 




Using these values of f\ and fa we can estimate the isospin-breaking correction 
to the energy shift of the 7r~p atom in the Is state. The strong energy-level 
shift t\ s of the -n~p atom is given by the model-independent formula [14]: ei s = 
e is° + e \ s LO = e L°(l ± fie) where the leading order (LO) or isospin-symmetric 
contribution is ef s ° and the next-to-leading order (NLO) or isospin-breaking 
contribution is e^ s LO . The quantity eff is expressed with the help of the well- 
known Deser formula [18] in terms of the S'-wave ttN scattering lengths with 
ef s ° = — 2a 3 fj J ^A s tr and A s tr = (2ai +03 )/3. The reduced mass of the n~p atom 
is denoted by fj, c = m p M n + /(m p + M 7r +) and A str = (88.4 ± 1.9) x 10 _3 M~+ is 
the strong (isospin-invariant) regular part of the n~p scattering amplitude at 
threshold [19] (for the definitions of these quantities see Ref. [14]). In ChPT 
the quantity 8 e , the ratio of NLO to LO corrections, is expressed in terms of 
the LECs ci, fa and fa 



= . *, MA t 8c i( M '+ " M *°) " e2 ( 4 /i + /*)] " 2«/i c (lna - l)A str . (32) 



The quantity C\ is the strong LEC from the ChPT Lagrangian [5,7] and F n = 
92.4 MeV is the physical value of the pion decay constant [14]. In Ref. [13] we 
obtained c x = —1.16 ±0.1 GeV -1 using the PCQM approach. Our prediction 
for C\ is close to the value C\ = —0.9mJ^ deduced from the irN partial wave 
analysis KA84 using Baryon Chiral Perturbation Theory [7] . Substituting the 
central values for our couplings fa = —19.5 MeV, fa = —8.7 MeV and c\ = 
—1.16 GeV -1 into Eq. (32), we get 5 t = — 2.8-10" 2 . Our estimate is comparable 
to a prediction based on a potential model for the irN scattering [19]: 



In conclusion, we give predictions for the 0(p 2 ) electromagnetic (e.m.) low- 
energy couplings (LECs) fa, fa and fa as originally set up in the ChPT effective 
Lagrangian. The magnitude of fa and its relation to fa and fa are obtained 
from an analysis of the nucleon e.m. mass shift and the leading radiative cor- 
rections to the irN scattering amplitude at threshold. Using our values for fa 
and fa we also predict the isospin-breaking correction to the strong energy 
shift of the n~p atom in the Is state. Latter prediction is extremely impor- 
tant for the ongoing experiment "Pionic Hydrogen" at PSI, which intends to 
measure the ground-state shift and width of pionic hydrogen (7r _ p-atom) at 
the 1% level [20]. 
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Fig. 1. Electromagnetic mass shift of the nucleon. 
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Fig. 2. Leading e /F radiative corrections to the ttN amplitude at threshold. 
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